Junction conditions at spacetime singularities 
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A classical model for the extension of singular spacetime geometries across their singularities 
is presented. The regularization introduced by this model is based on the following observation. 
Among the geometries that satisfy Einstein's field equations there is a class of geometries, with 
certain singularities, where the components of the metric density and their partial derivatives remain 
finite in the limit where the singularity is approached. Here we exploit this regular behavior of the 
metric density and elevate its status to that of a fundamental variable - from which the metric is 
constructed. We express Einstein's field equations as a set of equations for the metric density, and 
postulate junction conditions that the metric density satisfies at singularities. Using this model we 
extend certain geometries across their singularities. The following examples are discussed: radiation 
dominated Friedmann-Robertson- Walker Universe, Schwarzschild black hole, Reissner-Nordstrom 
black hole, and certain Kasner solutions. For all of the above mentioned examples we obtain a 
unique extension of the geometry beyond the singularity. 



I. INTRODUCTION 

General relativity (GR) is in excellent agreement with 
all experimental and observational tests of gravity. How- 
ever, when a GR solution develops a singularity it is 
sometimes impossible to use Einstein's field equations to 
extend the geometry beyond the hypersurface where the 
singularity resides. The standard lore is that this in- 
completeness should be resolved by a quantum gravity 
theory that may become important when the curvature 
length scale is of the order of Planck length. It is pos- 
sible that quantum-gravity phenomena can smooth out 
singularities and thereby resolve the nonextendibility of 
singular GR solutions. Alternatively, it is possible that 
the correct description of gravity in the vicinity of sin- 
gularities requires first an extension or a modification of 
classical GR (e.g. classical string theory introduces mod- 
ifications to GR [1]). In this case, quantization of gravity 
can take place only after our classical understanding of 
gravity has been modified Therefore, there is a good 
motivation to extend or modify GR so that it would not 
breakdown at singularities. Moreover, quantum theories 
can inherit singularities that appear in their correspond- 
ing classical theories (e.g. a charged point particle gives 
rise to an infinite electrostatic energy in classical electro- 
dynamics, a related ultraviolet divergent behavior also 
appears in quantum electrodynamics). Therefore, han- 
dling spacetime singularities at the classical level may 
provide a useful preliminary step towards a more regular 
quantum gravity analysis of singularities. Furthermore, 
this classical analysis could shed light on the extension 
of spacetime geometries beyond their singularities. 

In this manuscript we introduce a classical model for 
gravity that retains its predictive power for a class of 
spacetime singularities. In the construction of the model 
we make sure that the following requirements are sat- 



isfied: First, to ensure agreement with observational 
constraints, we demand that except for singularities the 
model would coincide with GR. Second, we demand that 
the model would be able to provide predictions for the 
extension of certain spacetime geometries beyond their 
singularities. 

Below we use our model to extend certain GR solutions 
across their singularities. The following examples are 
studied: Friedmann-Robertson- Walker Universe (FRW) 
, Schwarzschild black hole, charged Reissner-Nordstrom 
black hole (RN), and Kasner Universe. For other clas- 
sical approaches to singularities see Penrose [|| and Tod 

The regularization that we introduce is based on the 
following observation. There is a class of spacetime 
singularities where the components of the metric den- 
sity (in certain coordinates) attain finite values in the 
limit where the singularity is approached. As an ex- 
ample let us consider a Schwarzschild black hole. In 
the Schwarzschild coordinates the metric reads g^ h = 
diag(— b/r, r/b, r 2 , r 2 sin 2 9), where b — r — 2M, and 
M denotes the mass of the black hole. In the limit 
r — > 0, the Kretschmann scalar R a ^ lS R a p^s diverges 
thereby indicating that there is a spacetime singularity 
at r = 0. Here we focus attention on the contravari- 
ant metric density Q^ ch which is equal to the product 

[— det(g^ h )] 1 / 2 g^h- I n the Schwarzschild coordinates, 



it reads Q^ ch = diag[ 



1 sin 6, rb sin 8, sin 8, sin ( 
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Notice that in the limit of interest the components of 
dsch an d their partial derivatives remain finite, though 
the complete matrix Q§ ch becomes degenerate. In this pa- 
per we exploit this non-divergent behavior of the metric 
density and elevate its status to that of a fundamental 
variable, while the metric becomes a constructed vari- 
able, which is constructed from the metric density. We 
express Einstein's field equations as a set of equations for 
the metric density and postulate junction conditions at 
the singularities. We use this framework to extend cer- 
tain geometries across their singularities. At the singu- 
larities we find that the components of the metric density 
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are continuous and differentiable even though the geome- 
try remains ill-defined. In fact our model does not remove 
any physical singularity. On the contrary, as in GR we 
find that the Schwarzschild geometry becomes singular as 
r — > 0. However, in our framework the spacetime geom- 
etry becomes a constructed entity, which is constructed 
from a non-divergent fundamental quantity - the met- 
ric density. Since the metric is no longer a fundamental 
dynamical variable, its singular behavior at a spacetime 
singularity does not obstruct the continuation of the so- 
lution beyond the singularity. 

This paper is organized as follows: in Sec. [11] the clas- 
sical model is presented, in Sec. IHII we discuss few exam- 
ples in detail, and Sec. [TV] provides conclusions. 

II. THE CLASSICAL MODEL 

Let us begin by defining the metric density to be 
a four-dimensional symmetric contravariant tensor den- 
sity Q of weight -1. In the standard GR formulation the 
metric is represented by a non-degenerate matrix. Recall, 
however, that the metric density may become degenerate 
in the limit where a spacetime singularity is approached. 
For example, the Schwarzschild metric density Q<f ch be- 
comes a degenerate matrix in the limit r — > 0. Exploiting 
the non divergent behavior of g^ h , our first step is to ex- 
tend the solution in a continuous manner, and allow 
the metric density to attain its degenerate configuration 
in the limit. More broadly, a similar non-diverging be- 
havior of the metric density also appears in other GR 
solutions; for example, in appropriate coordinates, the 
components of the metric densities of RN, Kerr, Kerr- 
Newmann, FRW, and Kasner solution, all remain finite 
at their corresponding singularities (see discussion be- 
low). This means that within the space of GR solutions 
there is a class of singular solutions that approach fi- 
nite but degenerate configurations at their singularities. 
Our strategy is to include these degenerate metric density 
configurations in our formulation, and thereby construct 
a continuous extension of the space of GR solutions. 

Notice that it is possible to introduce a coordinate 
transformation such that in the new coordinates the met- 
ric density Qcf ch would diverge at the singularity. This 
divergent behavior is merely a coordinate singularity. In 
analogy with GR, we shall restrict attention to a class 
of nonsingular coordinates, where the metric density and 
its derivatives (see more details below) do not diverge. 

From the metric density we construct the contravariant 
metric through 

s-^H)- 1 ^. (i) 

Here 

= det(fl^) . (2) 

Except for being the quantity from which the metric is 
derived, the metric density also has its own independent 



physical meaning. In analogy with the metric, one can 
employ the metric density to construct scalars. These 
scalars are normally associated with a four-dimensional 
domain D. For example, the four volume of a domain D 
is defined by 

V= [ y/^Q^X. (3) 
JD 

The inner products between covector fields, and the inner 
product between vector fields are defined to be 

(a a ,bfj) = / a a b P Q aP d 4 x, (4) 

J D 

{u a y)= f u a v^ afS d 4 x. (5) 

J D 

Here the covariant metric density Q a p is defined to be the 
adjoint matrix (also called the adjugate matrix) of g Q/3 

0Q/3 = adj( ^), (6) 

This means that Q a p is equal to the transpose of the co- 
factor matrix of g Q °. From definition ([6]) it follows that 
Q a p transforms as a covariant tensor density of weight — 1 
under a general coordinate transformation [6|], and fur- 
thermore Eq. ^ guarantees that g Q ^ remains well de- 
fined even if Q a " becomes degenerate. For example, in the 
Schwarzschild singularity both q^q 1 and $f ch as well as 
their partial derivatives (up to all orders) , all have a well 
defined limit 1 as r — >■ 0. By contrast, the singular behav- 
ior of the Krctschmann scalar at the Schwarzschild sin- 
gularity implies that the metric must exhibit some kind 
of singular behavior at this limit 2 . 

Notice that " observables" (|3l4l5p that are constructed 
directly from the metric-density remain well defined even 
if domain D contains a singularity where takes the 
form of a nondiverging degenerate matrix. By contrast, 
a curve in D that passes through such singularity may 
have an ill defined length. From Eq. © we find that 

0°*W = fl^- (7) 

In analogy with the metric, one can use the metric 
density to construct evolution equations for other fields. 
Equations of motion that depend directly on the metric 
density may retain their predictive power at singularities 
where the metric density becomes degenerate. For ex- 
ample, in a fixed RN background the scalar wave equa- 
tion (Q a/3 (f>^a),p = determines the transmission of scalar 
waves through the RN singularity @- 



1 We ignore coordinate singularities, e.g. at 6 = 0, since these 
singularities can be removed by a coordinate transformation. 

2 More precisely, since the diverging Kretschmann scalar is equal to 
a sum of products of terms that depend on the metric (including 
the metric, its inverse, and their derivatives up to the second) 
it follows that in all possible coordinates at least one of these 
metric dependent terms must diverge at the singularity. 
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Next we define the covariant metric to be 

Sa/3 = -{-QY 1/2 Q a f3- (8) 

It follows from Eqs. (|1I2I7|8|) that for g^Owe have 




99af} 



a/3 _ 



(9) 
(10) 



Here g = det(g a p). 

We now construct the equations of motion for the met- 
ric density. For this purpose we cast Einstein's field equa- 
tions in a densitized format. Using the standard Landau- 
Lifshitz formulation H and units where G — c = 1 we 
obtain 



(0 



), a f3 
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Here T MI/ 
n vt<~ total 



(-0)T" V , where T^ v 
denotes the energy-momentum tensor depending on the 
matter fields and the metric density via Eqs. (|1I8[) . and 



'■ll 



1 



P 0Aq0 , p )/3 



0Aa0^0 i/3 F , p 



+ §0~ WV Q - 0"V Q )(20^fl CTT - P .0^)0^ A 0" CT C 

Here = (—0)tj^ L , where t^ L denotes the Landau- 
Lifshitz pseudo-tensor |9j. For g ^ we may substitute 
Eqs. ([TO]) into the equations of motion f|TT]) and recover 
the standard Einstein's field equations depending on the 
metric. In this case, Eq. (ITTI) implies that , satisfies 
the Landau-Lifshitz energy-momentum conservation law 



the components of the metric density would be analyt- 
ical functions of the coordinates. In fact requiring ana- 
lyticity of the metric is a standard method to continue 
the Kerr geometry through its ring singularity (see e.g. 
10]). However, in some respect demanding analyticity of 
the metric density is a too strong requirements. An ana- 
lytical continuation completely determines the extension 
of the solution, so it is no longer necessary to solve the 
equations of motion. This however contradicts our goal of 
keeping the standard physical picture where the solution 
is evolved via differential equations, and the junction con- 
ditions merely determines how the solution is extended 
across the hypersurface with the singularity. 

A weaker condition than analyticity is smoothness (i.e. 
C°°). Indeed an assumption of a smooth metric is often 
found in standard GR theorems. For example, in a rigor- 
ous initial value formulation one assumes smooth initial 
data [ll|. In analogy with this standard GR assump- 
tion, we postulate that the components of metric density 
be smooth functions of the coordinates. More precisely, 
we focus attention to metric density configurations that 
satisfy Eq. (fTTj) . not necessarily in vacuum, where g Q/3 
is smooth away from a genuine singularity of the geome- 
try, and its components and all of their partial derivatives 
have a well defined limit as the singularity is approached. 
We assume that g may be zero at most on a hypersur- 
face, but not in an open set (i.e. not in the bulk), and 
demand that at a singular hypersurface, where g = 0, the 
components of remain smooth. We then use these 
conditions together with Eq. (|TT|) to extend the solution 
across the singularity. Let us now discuss few examples 
in detail. 



III. EXAMPLES 
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Let us consider a vacuum singularity, e.g. the 
Schwarzschild singularity, where g Ml/ is finite but degen- 
erate. Notice that some of the terms in have nega- 
tive powers of g. These terms become ambiguous when 
the metric density becomes degenerate. This does not 
mean that t^ L must diverge at singularity, e.g. in the 
Schwarzschild coordinates v^ L has a well defined limit as 
the singularity is approached. Still the ambiguity in t^ L 
may give rise to difficulties in the extension of the solu- 
tion across the singularity by using Eq. (fTT|) alone. To 
overcome these difficulties we shall postulate new junc- 
tion conditions at the singularity. 

Our goal is to find junction conditions that can be 
combined with the equations of motion (|11[) in a man- 
ner that would provide a unique extension of the solution 
across the hypersurface where the singularity resides. We 
would like to exploit the fact that the components of the 
metric density and their derivatives may remain finite as 
the singularity is approached, and so it is possible to de- 
mand that these components would have some degree of 
smoothness at the singularity. It is tempting to demand 



A. Friedmann- Robertson- Walker 

First we consider a flat FRW Universe filled with a 
perfect fluid. In the standard GR formulation this solu- 
tion is nonextendible since its geometry becomes singular 
at the big bang. By virue of the underlying symmetry 
we express the metric density as g a/3 = A(rj)rj al3 , where 
T]^ denotes the Minkowski metric, and the scale factor 
is given by a = V~A, where A > 0. As in GR we fix 
the signature of the metric to be, say (— , +, +, +). The 
matter source term reads = dia,g(p,p,p,p), where 
p = A~ 3 p and p = A~ 3 p are the proper pressure and 
enregy density, respectively. For A ^ we substitute the 
expressions for T Ml/ and g M " into Eq. (fTTI) and obtain 



3i 2 = 327T/5 

3i 2 - AAA = 327if . 



(13) 
(14) 



Here an overdot denotes differentiation with respect to 
conformal time r\. We assume an equation of state of 
the form p = wp, where w is a constant. Substituting 
this relation into Eqs. (I13I14P gives p cx A 3( - 1 ^ w ' l/2 . We 
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substitute p(A) into Eq. (THfl) and solve for A. Recall 
that we seek a solution where Q a ^ and all of its par- 
tial derivatives have a well definied limit as singularity 
is approached. For the paricular coordinates in use, this 
condition is satsified for a radiation dominated Universe 
charterized by w = 1/3. Demanding that g Q/3 be smooth 
at the singularity gives 

A oaf. (15) 

This solution describes a bounce at the big bang. 

Notice that the matter source term is continu- 
ous at the singularity, but the complete source term 
^ total nas an ambiguity at A = 0, since contains 
the combination A/A, which is ambiguous for A = 0. 
Imposing continuity, we define t^(0) to be the limit 
of t^(rj) as rj — > 0. This gives, for all values of 77, 
167rt£i = -i 2 diag(3/2,7/2,7/2,7/2). Substituting this 
expression into the right hand side of Eq. ([TTj) shows that 
solution (TTSl) satisfies the equation of motion (fTTj) at the 
singularity. Moreover, evaluating the divergence of 2^ j 
shows that the energy-momentum conservation law (fTS)) 
is satisfied at the singularity. By contrast, it is hard to 
make sense of the covariant conservation law V Q T^" = 
at the big bang, since both V„ and T^ 1 diverge at the 
singularity. 

The above solution shows that the FRW singularity is 
described by a field configuration where both Q a @ and 
Q a /3 are smooth, but the metric as defined in Eqs. (jll8[) 
can not be constructed. This suggests the interpretation 
that a degenerate metric density describes a pre-metric 
configuration which is more primitive than a geometry. 

B. Black holes 

Next we study a spherically symmetric charged black 
hole characterized by a mass M and a charge Q. Sub- 
stituting the RN solution into Eq. (fit))) we find that for 
r > the non-zero components of the metric density are 
given by 

8 - = ^iw' A=sW (16) 

Srn = [Q 2 + r(r - 2M)] sin 6 , gf N = (sinOy 1 , 

and s/—0rn = r 2 smO. We focus on the domain r < r_, 
where r± = M ± y^ 2 - Q 2 , and M > |Q|. Notice 
that in the limit r — > the RN metric density becomes a 
non-divergent degenerate matrix. The r = surface cor- 
responds to a genuine singularity since the scalar R I1V R} 1V 
diverges in the limit. We now extend the solution across 
this singularity to the domain r < 0. For this purpose, 
we first have to solve Maxwell equations in vacuum. To 
overcome the singularity in the electromagnetic field we 
cast Maxwell equations in a densitized format, and ob- 
tain 

r^ = o. (17) 




FIG. 1: Penrose- Carter diagram of the extended RN space- 
time. 

Here S ra ' 3 = \J— QF a P , where F a/3 is the electromagnetic 
field tensor. For r > the non zero components of the 
RN solution 3#jy read 

3% N = -r£ N = Qsmd. (18) 

Notice that Eq. (|17J1 has no singularity at r = 0, and 
solution (|18p satisfies the Maxwell equations (|T7|) in the 
entire domain —00 < r < r_. We now solve Eq. (1111) in 
r < 0. Exploiting the underlying symmetry we substitute 
Q af3 = diag[-C 4 (r)sin6'/ J D(r),D(r)sin6',sin6',(sine')- 1 ] 
into Eq. (fTI]) . where we introduced the combination 
C 4 /D in g" to slightly simplify the equations. Demand- 
ing that g rr , g'* be smooth at r = 0, and using Eqs. (fTo) 
determines C(r) and D(r) uniquely. The complete so- 
lution for —00 < r < r_ is given by Eq. (|16j) . This 
expression for the extended geometry is also obtained 
for the special case of the Schwarzschild solution where 
Q = 0. Figure 1 shows the Penrose-Carter diagram of the 
RN extension for Q ^ 0. We should mention that this 
extension had been previously conjectured in Ref. Q. 
Notice that the components of this extended RN metric 
density are analytical in the domain of interest, and in 
particular they are well defined at the r — singular- 
ity. However, the metric as defined in Eqs. (|1I8[) can not 
be constructed at the singularity. Similar to the FRW 
example, we define 1^ tal (0) to be the limit of 3^ a jfV) 
as r — > 0. With this definition, the equation of motion 
(TlTj) and the conservation law (TT2"|) are satisfied at the 
singularity. 

The geometry of a rotating Kerr (and Kerr-Newmann) 
black hole can be analytically extended through its ring 
singularity using the standard GR formalism (see e.g. 
[lOj]), and therefore a regularization of the singularity 
is not essential in this case. Nevertheless, it interest- 
ing to note that in the Boyer-Lindquist coordinates the 
components of the metric density of the Kerr (and Kerr- 
Newmann) solution are analytical in the neighborhood of 
the singularity. Here again the metric density takes the 
form of a degenerate non-diverging matrix at the singu- 
larity. 
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C. Kasner 

The Kasner solution is a vacuum solution which is flat 
and homogenous but anisotropic. In the Kasner coordi- 
nates (t, x, y, z) the metric density of the Kasner solution 
is given by 

gf as = diag(-t, t 1 - 2 " , t 1 - 2 ^ , i 1 " 2 ^ ) . (19) 

Here the three parameters pi, satisfy 

Vi+P2+P3=p\+Pl+Pl = l- (20) 

The Kasner geometry is regular for t > but has a sin- 
gularity at t = 0. To be able to use our model and extend 
the geometry to the domain t < 0, we should check that 
the components of the metric density, and their partial 
derivatives, have a well defined limit as the singularity 
is approached. Here however there is a difficulty. Notice 
that the constraints (|20[) imply that at least one of the pi 
parameters is greater than 1/2, and so at least one of the 
components of Q^f as must diverge at the singularity. To 
resolve this difficulty we should seek a coordinate trans- 
formation that removes the singularity from Qjf as - Below 
we provide such a transformation for the case where the 
Pi parameters are given by three different rational num- 
bers. 

Let us consider first the domain t > and introduce 
the coordinate transformation t = eri s , where s is a posi- 
tive integer whose value is specified below, and e = 1. In 
the new coordinates (77, x, y, z) the metric density reads 

&ial = diag(-s- V V 1 , V 2 , V 3 ) • (21) 

Here qi = 2s(l — pi) — 1. Let us denote the combinations 
(1 — pi) with rrii/rii, where mi,m are positive integers. 
Setting s — n\nin^ implies that the parameters 91, 92, Q3 
are given by positive odd integers. Notice that the com- 
ponents of the Kasner metric density 0^f s and their par- 
tial derivatives, both have a well defined limit as r\ — > 
as desired. 

We are now ready to use our model and extend the 
solution to 77 < 0. Notice that solution ([21]) satisfies Eq. 
(ITU in 77 < 0, though here the signature is (+,—,—,—). 
Imposing our junction conditions by demanding that the 
components of the metric density be smooth at the r\ = 



singularity, ensures that the parameters s and q\ , qi , qz 
have the same values for 77 < and 77 > 0. This uniquely 
determines the continuation of the geometry across the 
singularity. It is now possible to transform back to the 
original coordinates (For the case where s is an even num- 
ber, we set e = —1 for 77 < 0, so that negative values of 
77 correspond to negative values of t) . 



IV. CONCLUSIONS 

Singularities mark the breakdown of the laws of gen- 
eral relativity. There is a widely accepted view that a 
viable quantum theory of gravity should have a mech- 
anism that removes singularities, for example by some 
quantum-gravity phenomena that smooth them out. In 
this paper we have presented a classical model that illus- 
trates an alternative point of view. We showed that it is 
possible to keep the singularities in the solution, so that 
they become a legitimate part of the physical model. 

The construction of our model was based on the obser- 
vation that the metric density encodes the same amount 
of information as the metric, but unlike the metric, its 
components and their partial derivatives may remain fi- 
nite at a class of spacetime singularities. Our strategy 
was to exploit this property, and replace the metric with 
the metric density as the fundamental variable of gravity. 
By retaining the dynamical evolution of GR, and supple- 
menting it with junction conditions we showed that it 
is possible to extend certain singular geometries across 
their singularities. 

It would be interesting to see if our model could be 
used to extend singular geometries in d dimensions. Cu- 
riously, however, not all the observables can be adjusted 
to accommodate a d dimensional metric density. In par- 
ticular notice that Q a p transforms as a tensor density of 
weight — (d— 3), and so the inner product (m q ,v' 3 ) de- 
fined by Eq. ([5]) is a scalar only if d = 4. 
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